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Anomalous phonon red shift in K-doped BaFe2As2 Iron-pnictides 
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The effect of K-, Co- and P-doping on the lattice dynamics in the BaFe 2 As 2 system is studied by 
infrared spectroscopy. We focus on the phonon at ~ 253 cm~^, the highest energy in-plane infrared- 
active Fe-As mode in BaFe 2 As 2 . Our studies show that the Co- and P-doping lead to a blue shift 
of this phonon in frequency, which can be simply interpreted by the change of lattice parameters 
induced by doping. In sharp contrast, an unusual red shift of the same mode was observed in the 
K-doped compound, at odds with the above explanation. This anomalous behavior in K-doped 
BaFe 2 As 2 is more likely associated with the coupling between lattice vibrations and other channels, 
such as charge or spin. This coupling scenario is also supported by the asymmetric line shape and 
intensity growth of the phonon in the K-doped compound. 

PACS numbers: 74.25.Gz, 78.30.-j, 74.25.Kc 


In spite of extensive studies on high-Tc superconduc¬ 
tivity in iron pnictides since its discovery^i the question 
of what plays an important role in the paring mechanism 
of this class of superconductors remains enigmatic. The¬ 
oretical calculations have demonstrated that, unlike the 
traditional BCS superconductors, an electron-phonon in¬ 
teraction is not sufficient to account for such a high 
in iron pnictidesi^ Therefore, the s± paring with a sign 
reversal in the gap function, that is mediated by spin 
fluctuations, was proposed by Mazin et al.^ However, 
the above scenario seems far from being thoroughly es¬ 
tablished. The s± paring state is expected to be sensitive 
to impurity scattering/ Indeed, on one hand in- and out- 
of-plane dopings affect T), and the unpaired quasiparticle 
density in the BaFe 2 As 2 family^ but on the other hand, 
superconductivity is robust against impurities in 1111 
materialsFurthermore, a large iron isotope effect has 
been reported in SmFeAs 0 i_a;Fe 2 , and Bai_2,K3;Fe2As2r 
indicating that electron-phonon coupling plays some role 
in the paring mechanism. By taking these facts into 
account, Kontani et al. proposed that electron-phonon 
coupling arising from the Fe-ion oscillation can induce or¬ 
bital fluctuations, mediating the s++ paring without sign 
reversal.— A recent quantitative convergent beam electron 
diffraction study by Ma et aL^ has revealed strong cou¬ 
pling between Fe orbital fluctuations and anion-dipole 
polarizations in Ba(Fei_a;Coa;) 2 As 2 . They suggest that 
a full understanding of the paring mechanism in iron- 
pnictides can only be reached by considering the charge, 
spin, orbital, lattice and anion polarization all together 


in a consistent theory. 

Infrared spectroscopy is a standard tool to investigate 
lattice dynamics, providing information on the coupling 
between lattice vibrations and electrons or spinsj^i// Al¬ 
though optical investigations into the lattice dynamics of 
iron pnictides have been conducted by many groups 
a comparison study of different iron pnictides, in partic¬ 
ular different substitution types in the same family, has 
never been performed. 

We fill this gap by comparing the behavior of the 
^253 cm“^ in-plane infrared-active Fe-As mode in 
BaFe 2 As 2 (parent compound of the Bal22 family) and 
three different doping types (K, Co and P doping). We 
observe a blue shift of this mode in the Co- and P-doped 
BaFe 2 As 2 , as well as an anomalous red shift of the same 
mode in the K-doped compound. The latter can not be 
explained by the change of lattice parameters induced by 
doping. A close inspection of the phonon line shape and 
intensity leads us to the conclusion that the coupling be¬ 
tween lattice vibrations and other channels is stronger in 
K-doped BaFe 2 As 2 and thus responsible for the anoma¬ 
lous red shift of the phonon. 

High quality single crystals of BaFe 2 As 2 (BFA, 

~ 138 K), Bao.6Ko.4Fe2As2 (K40, Tc ~ 39 K), and 
Ba(Feo, 92 Coo,o 8)2 As 2 (Co08, Tc ~ 23 K) were grown with 
a self-flux method /®/® while the BaFe 2 (Aso, 85 Po.i 5)2 
(P15, Tjv — 90 K) single crystals were grown without 
fluxj^ The a6-plane reflectivity R{ui) was measured at 
a near-normal angle of incidence on Bruker IFS113v 
and IFS66v spectrometers. An in situ gold overfilling 
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Figure 1. (color online) The reflectivity in the far infrared 
region for BFA, K40, Co08 and P15 at 300 K (a) and 150 K 
(b). The red arrows indicate the afe-plane infrared-active Ba 
mode at 94 cm“^ and Fe-As mode at 253 cm“^ in BFA. 


technique^! was used to obtain the absolute reflectivity of 
the samples. Data from 30 to 15 000 cm“^ were collected 
at different temperatures on freshly cleaved surfaces for 
each sample, and then we extended the reflectivity to 
40 000 cm“^ at room temperature with an AvaSpec-2048 
X 14 optical fiber spectrometer. 

Figure [IJa) shows the far-infrared reflectivity at room 
temperature (300 K) for all 4 compounds. Metallic be¬ 
havior can be realized in all these materials by their rel¬ 
atively high reflectivity that approaches unity at zero 
frequency. In addition, two sharp features (indicated 
by the arrows), representing the symmetry-allowed ab 
plane infrared-active modes, were observed at ~ 

94 cm“^ and 253 cm“^ in BFA, consistent with previous 
works j^^d^ The 253 cm“^ mode is also clearly observed 
on the room temperature R(oj) in the doped BFA com¬ 
pounds. Figure[TKb) displays R{uj) for the same materials 
at 150 K, where similar features are revealed. 

Here we would like to point out that BFA exhibits 
structural and magnetic phase transitions at 140 K— 
with a consequent renormalization of the infrared phonon 
spectra^^— whereas such transitions are either sup¬ 
pressed or absent in the doped compounds, meaning that 
below 140 K, these compounds are in different phases. In 
order to avoid effects related to these phase transitions, 
the temperature window in our study is constrained be¬ 
tween 150 K and 300 K. 

In the following, we concentrate on the 253 cm“^ 
mode, which involves the displacements of Fe and As 
atoms j^^d^ To investigate the doping effect on this mode 


Figure 2. (color online) The low frequency optical con¬ 
ductivity cri(a;) for BFA, K40, Co08 and P15 at 300 K (a) 
and 150 K (b). Panel (c) and (d) show the enlarged view 
of Panel (a) and (b), respectively, focusing on the phonon at 
253 cni“^(with offset). The thin black lines through the data 
are the Lorentzian fitting results. The vertical dashed lines 
indicate the phonon frequency for BFA. 


in a straightforward way, we calculated the optical con¬ 
ductivity via Kramers-Kronig analysis of the reflectiv¬ 
ity. At low frequency, we employed a Hagen-Rubens 
{R = 1 — extrapolation. Above 40 000 cm~^ (the 

highest measured frequency), we utilized a constant re¬ 
flectivity up to 12.5 eV, followed by a free-electron (w^) 
response. 

FigurellKa) and Fig.[^b) show the real part of the op¬ 
tical conductivity ai(uj) in the far-infrared region for all 
compounds at 300 K and 150 K, respectively. The low- 
frequency CTi (w) exhibits a prominent Drude-like metallic 
behavior for all materials at both 300 K and 150 K, con¬ 
sistent with the reflectivity analysis. 

The optical conductivity can be conveniently parame¬ 
terized by a Drude-Lorentz model: 
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where Zq is the vacuum impedance. The first term de¬ 
scribes a sum of free-carrier Drude responses, each char¬ 
acterized by a plasma frequency flpj and a scattering 
rate 1 /r^. The second term is a sum of Lorentz oscilla¬ 
tors, each having a resonance frequency wo.fc) a line width 
7 fc and an oscillator strength Hfe. The optical response of 
FeSCs can be modeled reasonably well by the superpo¬ 
sition of two Drude components and a series of Lorentz 
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Figure 3. (color online) The resonance frequency cjq of the 
Fe-As mode as a function of temperature for all the com¬ 
pounds. The black dashed lines represent a quadratic tem¬ 
perature dependence. 


terms over a wide frequency range, which has been dis¬ 
cussed in detail in previous works<^“— 

In addition to the gross features, the 253 cm“^ mode 
manifests itself as a sharp peak in the optical conduc¬ 
tivity. Figure [5Kc) highlights the region around the 
253 cm“^ phonon for all compounds. The dashed line 
denotes the phonon peak position of BFA. It can be im¬ 
mediately noticed that the phonon in the K-doped Bal22 
compound shifts to lower frequency (red shift or soften¬ 
ing) , while P- and Co-doping lead to a shift of this mode 
to higher frequency (blue shift or hardening). Exactly 
the same behavior is observed at 150 K, as shown in 
Fig.IKd). 

In order to quantitatively analyze the behavior of the 
phonon upon doping, we fit it to a Lorentz oscillator with 
a linear background in a narrow frequency range centered 
at the phonon resonance frequency for all the materials at 
all measured temperatures. The fitting results at 300 K 
and 150 K are shown as thin solid lines through the cor¬ 
responding data in Fig. [5](c) and Fig. [5Kd), respectively. 
The fitting parameters are summarized in Table ID 

Figure [3] shows the phonon resonance frequency ojq de¬ 
termined from the fit for all the compounds at 7 measured 
temperatures between 150 K and 300 K. For each mate¬ 
rial, ujQ increases upon cooling, following a quadratic T 
dependence (indicated by the black dashed lines), which 
is expected in the absence of structural or magnetic tran¬ 
sitions. The T dependence of this phonon is in good 
agreement with previous works For all the measured 
temperatures, with respect to BFA (red solid squares), 
cjQ is smaller in K-doped compounds (blue solid circles), 
but larger in the Co- (green solid diamonds) and P-doped 
(pink solid triangles) materials. 

We now trace the origin of the phonon frequency shift 
induced by different doping types. Doping usually has 
three main effects on materials: (i) changing the carrier 
concentration by adding electrons (electron doping) or 
holes (hole doping); (ii) introducing disorder; (hi) appli- 


Table I. The vibrational parameters for oscillator fits to the 
infrared-active Fe-As mode observed in different doped com¬ 
pounds at 300 K and 150 K, where ujo, D, and 7 are the os¬ 
cillator frequency, strength, and line width, respectively. All 
units are in cm“^. 


Doping 


300 K 



150 K 


LJQ 

Q 

7 

Wo 

Q 

7 

BFA 

253.10 

211.85 

3.77 

258.00 

217.13 

3.43 

K40 

251.30 

321.50 

5.51 

255.88 

322.74 

4.51 

C 0 O 8 

254.94 

208.26 

6.31 

259.75 

208.86 

5.93 

P15 

257.34 

175.10 

5.92 

261.78 

181.78 

5.88 


cation of chemical pressure that can change the lattice 
parameters. Increasing the carrier concentration will in¬ 
crease the electronic background, potentially leading to 
a growth in the screening effect. However, the screening 
effect only affects the intensity of the phonon, but has 
no influence on the frequency;^ In addition, the screen¬ 
ing effect on the 253 cm“^ mode in the Bal22 sys¬ 
tem has been demonstrated to be very weak or absent .f- 
Therefore, we can safely rule out this possibility. Dis¬ 
order in the Fe-As layers may reduce the intensity of 
this mode, but virtually has no effect on the phonon 
frequency either. Finally, the change of lattice param¬ 
eters due to chemical pressure may play an important 
role in shifting the phonon frequency. Generally, the fre¬ 
quency shift of a phonon is closely related to the change 
of the bond length I of the associated atoms. The phonon 
shifts toward higher frequency if I shrinks. The expected 
frequency shift can be estimated via a simple formula 
loo/ujq = (I 'By comparison with BFA, the Fe- 
As bond length I is 0.3% shorter in 0.5% 

shorter in Co08^ and 0.8% shorter in PI5^ Based on 
these parameters, the doping-induced frequency shift of 
the 253 cm“^ mode in BFA can be easily estimated in all 
three compounds. The open diamonds in Fig.^Ka) denote 
the calculated phonon frequency by considering the Fe-As 
bond length for each material, which we compare to the 
values determined from the experiment (solid squares). 
Note that the measured phonon frequency agrees very 
well with the calculation in both C 0 O 8 and P15, suggest¬ 
ing that the phonon frequency shift is dominated by the 
Fe-As bond length. In sharp contrast, a striking differ¬ 
ence between the calculated phonon frequency and the 
experimental result is found in K40. K doping induces 
a contraction of the Fe-As bond and, hence, one expects 
a blue shift of the corresponding phonon. However, a 
prominent red shift of the phonon frequency is observed 
from the experiment, which can not be explained by the 
change of the Fe-As bond length alone. 

This brings us to the question of what causes the 
anomalous red shift of the phonon in K40. One possi¬ 
bility is that the phonon is coupled to other channels, 
such as charge or spin. In the quasi-two-dimensional 
quantum spin system Sri_ 2 :Baa;Cu 2 (B 03 ) 2 , strong spin- 
phonon coupling leads to an additional 3% softening of 
the related phononj ^^i^^ The pronounced gate-induced 
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phonon softening in bilayer graphene is attributed to 
the coupling of the phonon to electronic transitions 
In addition, electron-phonon and spin-phonon coupling 
are widely observed in FeSCsi^^— Theoretical calcula¬ 
tions also point to spin-phonon coupling as the origin 
of the phonon softening in FeSCs .^^do This makes the 
electron- or spin-phonon coupling scenario more favor¬ 
able to accounting for the anomalous phonon red shift in 
the K-doped compounds. 

Further evidence for electron- or spin-phonon coupling 
can be revealed from the line shape of the phonon, since 
coupling of the lattice vibrations to charge or spin excita¬ 
tions often results in an asymmetric line shapei ^^d^ One 
may already notice that, as shown in Fig. [2jd), while a 
Lorentz oscillator yields a reasonably good description to 
the phonon line shape in BFA, Co08 and P15, the Lorentz 
fitting result for K40 is relatively poor: the phonon in 
K40 exhibits a slightly asymmetric line shape. We can 
emphasize the asymmetry of phonon line shapes by sub¬ 
tracting the background. The inset of Fig.Hl^b) shows the 
phonon for all compounds, with the background removed 
from the total optical conductivity. The asymmetric line 
shape in K40 (blue solid circles) can be easily distin¬ 
guished from other three materials, suggesting stronger 
electron- or spin-phonon coupling in the K-doped com¬ 
pound. 

The asymmetry of the phonon can be quantified by 
fitting it to a Fano line shape^i 


CTl (w) 


2 ^ + 

Zo 7 ’ 


( 2 ) 


where ojq, 7 and 0 represent the resonance frequency, 
line width and strength of the phonon, respectively. The 
asymmetry of the Fano line shape is described by a 
dimensionless parameter 1/q^. As \/q^ increases, the 
asymmetry of the line shape intensifies, indicating a 
growth of the coupling strength. In the case 1/q^ = 0, 
the symmetric Lorentz line shape is fully recovered. The 
inset of Fig. HKb) depicts the fitting results at 150 K for 
each material. The parameter \/q^, determined from the 
Fano fit, is shown as green solid circles in Fig.|4Kb). While 
the Fano parameter {1/q^ ~ 0.0016) is vanishingly small 
in BFA, C 0 O 8 and P15, its value (1/g^ ~ 0.045) is almost 
30 times larger in K40. This implies that the asymme¬ 
try of the phonon line shape is much stronger in K40, 
and thus the coupling between the phonon and electrons 
or spins is significantly enhanced in the K-doped com¬ 
pounds. 

The coupling between lattice vibrations and electrons 
or spins could also cause an increase of the phonon 
intensityThe intensity of the phonon deter¬ 
mined from the Fano fit at 150 K, is traced out for all the 
compounds as blue solid triangles (error bars within the 
symbols) in Fig. |4Kb). Comparing to BFA, 17^ doubles in 
K40, but slightly decreases in C 0 O 8 and P15. The change 
of the phonon intensity (the area under the phonon line 
shape) can also be identified directly from the phonon 



Figure 4. (color online) (a) The black solid squares denote 
the resonance frequency ujq of the Fe-As mode at 150 K for 
BFA, K40, C 0 O 8 , and P15. The red open diamonds are the 
calculated frequency by considering the Fe-As bond length 
change only for each compound, (b) The green solid circles 
and blue solid triangles portray the Fano parameter 1/q^ and 
phonon intensity determined from the Fano fit at 150 K, 
respectively. Inset: the phonon line shape without the elec¬ 
tronic background for all the compounds. The solid lines 
through the data are the Fano fitting results. 


line shape as shown in Fig. HKc), Fig. [Kd) as well as the 
inset of Fig.HKb). The slight decrease of 17^ in C 0 O 8 and 
P15 may arise out of disorder effect induced by in-plane 
doping. The striking increase of the phonon intensity in 
K40 is supportive of coupling between the phonon and 
electrons or spins in the K-doped compounds. 

One should note that our observations do not conflict 
with the evidence for the electron- or spin-phonon cou¬ 
pling in BaFe 2 As 2 Co-doped BaFe 2 As 2 ^ or other iron- 
pnictideS)^^— because our studies are based on the com¬ 
parison between the parent compound BFA and doped 
materials. We are not saying that there is no coupling 
in BFA and Co- or P-doped materials. What we assert 
is that the coupling is stronger in the K-doped material 
when compared to the other compounds. The relatively 
stronger electron- or spin-phonon coupling might be an 
important factor for the relatively higher Tc in the K- 
doped BaFe 2 As 2 system. 

In summary, the effect of three different substitution 
types on the 253 cm“^ a 6 -plane infrared-active Fe-As 
phonon in BaFe 2 As 2 has been studied by infrared op¬ 
tical spectroscopy. This mode shifts to higher frequen¬ 
cies in Co- and P-doped BaFe 2 As 2 , in agreement with 
the expected frequency shift associated with lattice pa¬ 
rameter change. Intriguingly, the same mode exhibits 
an unusual red shift in the K-doped compound. We at- 
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tribute the anomalous red shift of the Fe-As mode in 
K-doped BaFe 2 As 2 to electron- or spin-phonon coupling. 
The asymmetric line shape and intensity growth of this 
mode in K-doped BaFe 2 As 2 also support our conclusion. 
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